Topological Insulating States in Ordinary Semiconductors 
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We propose that ordinary semiconductors with large spin-orbit coupling (SOC), such as GaAs, 
can host stable, robust, and tunable topological states in the presence of quantum confinement 
and superimposed potentials with hexagonal symmetry. We show that, unlike other systems, the 
electronic gaps can be as large as the electronic bandwidth and that chiral spin edge states only 
exist because of the SOC. 
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A topological insulator (TI) is a fascinating state 
of matter which presents unusual physical properties 
such as a quantum Spin Hall effect in two-dimensions 
(2D), spin-polarized chiral Dirac surface states in three- 
dimensions (3D), exotic magneto-electric effects, and Ma- 
jorana fermions in the presence of superconductivity [lj- 
3]. Although TIs have attracted a lot of attention, 
progress in this area has been hindered by the absence of 
experimental systems with robust electronic and struc- 
tural properties. Although these materials have been 
christened as "topologically protected", material issues 
associated the weak strength of the SOC, the small size of 
the gaps, and the strong disorder (of the order of the elec- 
tronic bandwidth) present in most of the proposed sys- 
tems, make the experimental realization of these amazing 
physical properties very difficult, if not unattainable. 

One of the first theoretical predictions for such TI 
states was made for graphene [J] in the early days of 
graphene research Nevertheless, carbon is a light 
element with weak intrinsic SOC (s» 10~ 2 cV). In addi- 
tion, in flat graphene, due to wavefunction orthogonal- 
ity, the SOC for the 7r-bands is even weaker (ss 10 -6 eV) 
and, therefore, essentially unobservable. While devia- 
tions from flat sp 2 , to out of plane sp 3 , bonds can lead to 
three fold enhancement of the SOC [6(, the atomic con- 
trol over these deformations is a major experimental chal- 
lenge. Shortly after its initial proposal, TI states where 
predicted to occur in HgTe quantum wells ml, in 3D bulk 
solids of binary compounds involving Bi [8], and half- 
Heusler ternary compounds Nevertheless, all these 
materials have complex crystal structure which is very 
sensitive to stoichiometry. Hence, the unavoidable pres- 
ence of defects, which are usually unitary scatterers and 
can act as donors/acceptors, has a strong effect in the 
electronic structure. One observes, for instance, large 
broadening of the spectral lines for surface states (short 
lifetimes) in angle resolved photo-emission (ARPES), and 
doping of the bulk crystal ultimately transforming the 
TIs into metals [Icj . 

In this paper, we propose an alternative way to pro- 
duce robust, structurally stable, and tunable TI states in 



ordinary semiconductors such as GaAs. The advantages 
of these materials are that they have significant SOC, 
they can be grown with extreme precision using molecu- 
lar beam epitaxy (MBE) with large electronic mobilities, 
they can be tailored into quantum wells with arbit rary 
thickness, and can be controlled by external gates 11 1. 
The robustness and flexibility of these systems can be 
the starting point for the creation of different kinds of 
TIs that cannot be obtained otherwise. 



In fact, it is known that hole-doped zinc-blend semi- 
conductors naturally have large SOC that originates from 
the atomic p 3 / 2 ~ P1/2 nne structure splitting. Since the 
hole energy is always smaller than the splitting, the SOC 
is larger than the Rashba and Dresselhaus couplings that 
are relatively weak and of the same order of magnitude 
for both electrons and holes. We demonstrate that the 
effective SOC in a semiconductor quantum well with a 
superimposed hexagonal superlattice can be controlled 
by the strength of the transverse confinement and the 
scale of the superlattice. Hence, the SOC gap can made 
comparable to the bandwidth or continuously switched to 
zero. Finally, we show that the SOC leads to the appear- 
ance of chiral spin edge modes in contrast with systems 
such as graphene where these edge states exist even in 
the absence of SOC jl2j. Thus, in our proposal, the sys- 
tem can be continuously tuned between the Dirac metal, 
topological insulator, and standard band insulator. 



As it is well known [ll| , in the 3D bulk of a semicon- 



ductor, the hole wave function in systems like GaAs orig- 
inates from the atomic p 3 / 2 orbital and thus, the hole has 
an angular momentum J = 3/2 (the so-called, hole spin 
S = 3/2). In the large wavelength approximation (the 
k • p approximation), the hole effective Hamiltonian is 
proportional to the second power of the hole momentum 
k. The only kinematic structures allowed by symmetries 
are k , {k • S^j , and T^i/^^k^ki/S^S-y, where T^^^^ is the 
4 th rank tensor built of unit vectors of the cubic lattice 
that is known to be parametrically small relative to the 
other terms [lH and will be disregarded in what follows. 
In this case, the effective Hamiltonian can be written as 
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(we use units such that h = 1): 



with states given by: 



ff 3D = ^-(7i + 2.57)-7^-^, (1) 
2m e m e 

where m e is the free electron mass and 71, 7 are the 
Luttinger-Kohn parameters 13J. Hereafter, we use the 
parameters for GaAs where 71 ~ 6.8, 7 « 2.9 Notice 
that 7 parametrizes the effective SOC and is compara- 
ble with 71 , which parametrizes the effective hole kinetic 
energy. 

The 3D semiconductor can be geometrically confined 
in one direction creating at 2D quantum well. For sim- 
plicity, we assume that the confinement is described by 
an infinite square well of width d and that only the lowest 
quantum state, |0), has to be taken into account. In this 
case we have (0|fe z |0) = 0, and (0|fcf|0) = ir 2 /d 2 . As a 
result, the in-plane momentum is small k x ,k y <C ir/d 
and the k 2 S 2 term in {1} enforces the spin quantiza- 
tion along the z-axis. The lowest energy state corre- 
sponds to S z = ±3/2 and the higher state corresponds 
to S z = ±1/2 giving rise to the heavy (S z = ±3/2) and 
light (S z = ±1/2) hole states According to © the 
value of the splitting between these states is given by: 

A = E 1/2 - E 3/2 = 47^^ . 
' ' 2m P 
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Here we have introduced an effective "spin" s = 1/2 de- 
gree of freedom describing the two | f) and | I) states. 
Interestingly, 7, which can be considered as a strength of 
the spin orbit interaction, see Eq. ([T]) , is cancelled out in 
the expression for a in Eq. ([3]). One can call it "ultra- 
relativistic" behaviour, the spin-orbit is so large that it 
does not appear explicitly in the answer. 

In order to generate the TI, a potential U(r) with 
hexagonal (triangular) symmetry and spacing L, as 
shown in Fig[T] is superimposed to the 2D electron gas 
[l5j |. The lattice translation vectors are L\ = (L,0) , 
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When the hole density is low only the heavy hole band 
is filled. The heavy hole 2D dispersion follows from (fTJ) 
and has several contributions. Firstly, there is a diagonal 
term due to the z-confinement contribution which is just 
the matrix element of (fTJ): 



P 

2m, 



(7i + 2.5 7 ) 



- P 

■ 7 



2m, 



■(3/2|5 2 |3/2) = £-( 7 i+7) (2) 



Here p — {k x ,k y ) is the in-plane momentum. There 
is also the 2 nd order perturbation theory contribution 
due to the -^-[k x k z {S x S z + S z S x ) + k y k z (S y S z + S z Sy)} 
term in (flj. This term generates virtual z-excitations. A 

straightforward calculation gives the following 2 nd order 

2 

contribution for GaAs: — l-Q-^T- Putting these results 
together with © we find that the in-plane mass of the 
heavy hole is given by: 



m* — m e j (71 + 7 — 1.6) « O.I2m e . 

For a soft parabolic z-confinement the mass can be some- 
what larger, m* ~ 0.15 — 0.2m e . 

The (p ■ S) 2 part of the Hamiltonian (JT]) leads to the 
heavy-light hole mixing. The mixing matrix elements 
are: 

(p, -l/2\H\p, 3/2) = + Wyf 



(p,l/2|ff|p,-3/2) 



y/37 
2m e 



(Px ~ ip y ) 2 , 



FIG. 1: 



L 2 = (■§-, 2^)- Hence, there are two independent recip- 
rocal lattice vectors in the Brillouin zone (see Fig[T]) : 

2tt 2tt 
Gi = — (3,73), G 2 = — (0,2V3), G 3 = G X -G 2 , 

The points K \ . K 2 , K3 are connected by vectors G;, 
and K[ are obtained from the Ki by reflection. In order 
to simplify notation, we will measure energy in units of 
the bandwidth: 



En = 



K 2 
2m* 



(4tt/3£) s 
2m* 



In the case of GaAs, assuming L — 20 nm, which can 
be obtained experimentally with standard lithographic 
techniques, we have En 13 meV. Notice, however, this 
energy scale can be easily controlled by tuning L (for 
L = 50 nm, E « 2 meV. 

Unlike the case of graphene where the starting point 
is a tight-binding description [5], our description starts 
from a nearly free electron description. We assume, for 
simplicity, a periodic potential with a single Fourier com- 
ponent: 

U(r) = 2W [cos(Gi • r) + cos(G, • r) + cos(G 3 • r)} , 

(4) 

where W gives the strength of the potential. This poten- 
tial has nonzero matrix elements only between states |fe) 
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and \k ± Gi) with matrix elements given by W . Diago- 
nalization of the Hamiltonian, 



H 



P 

2m* 



U{r) 



(5) 



gives the hole dispersion shown in Fig.[2ja) with the pres- 
ence of two Dirac points with linear dispersion. If the 
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FIG. 2: (a) The hole dispersion along a particular contour in 
the BZ, (b) Map of the total charge density (in units 1/L 2 ) at 
the chemical potential tuned to the Dirac point. Both figures 
correspond to W = W/Eo = 1 in the absence of SOC. 

chemical potential is tuned at the Dirac point the av- 
erage hole density is (n) = 8/(3i 2 ). The map of the 
charge density is shown in Fig. &b) . At L = 50nm the 
average density is 1.1 x 10 n cm~ 2 . Even when the po- 
tential is strong, namely W = W/Eq = 1 — 2, the disper- 
sion, Fig. HJa) , is rather close to the result obtained by 
perturbation theory. The charge density plot, Fig. EJb), 
is fully connected with empty spots at positions of the 
potential maximums. So, in clear contrast to graphene, 
a t W ^ 2 the system is much closer to the nearly free 
electron regime than to the tight-binding one. 

Perturbative analysis of the system in the nearly free 
electron regime, at small W — W/Eo, is straightfor- 
ward [lij . A hole state close to the Dirac point, g« 1, 
is described by degenerate perturbation theory as: 



i> q cx ci|l) +c 2 |2) +c 3 |3) 



= i(*: 3 +q)-r 



(6) 



In the basis of states \j) the Hamiltonian ([5]) is repre- 
sented by 3 x 3 matrices: 



2m* 



2m* 









{K 3 +q) 



E 





(K 2 +q) 2 

2m* 



2m* 

K 1 q 
K 2 ■ q 
K 3 ■ q 




Notice that the eigen-energies of the U-matrix are 
0,0, 3W. In order to project in the double degenerate 



subspace of U we define [151 ]: 




i 

V2 



\b) 




(7) 



Projecting the kinetic energy to this basis and shifting 
the zero energy level to Eq one finds: 



(a\H\a) -> 
(b\H\b) = 
(b\H\a) = 



(K 2 + K 3 ) 
2m* 
(4Jd + K 2 



-VQx , 
Q 



(K 3 



6m* 
K 2 )-q 



+vq x 



where 



2V3r, 



K 



= -vq v 



2m* 



2vr 
3Lm* 



(8) 



is the Fermi-Dirac velocity. Notice that the velocity is 
controlled by the lattice spacing. Hence, in the Pauli 
matrix (pseudo-spin) representation the effective low en- 
ergy Hamiltonian reads: 



H 



v (-cr z q x 



(9) 



One can perform the unitary transformation H — > 
T^HT, where T represents two subsequent 7r/2 rotations 
around x- and y-axes in the pseudo-spin space and trans- 
form the Hamiltonian to the conventional form of a 2D 
Dirac Hamiltonian: H —> vcr ■ q. However, in what fol- 
lows we will use ©, as it is slightly more convenient for 
the study of the edge states. 

In order to account for the SOC one needs to intro- 
duce the light hole admixture in wave functions (|3]). This 
admixture leads to the following SOC correction to the 
matrix element of the potential 



S((P2\U\ Pl )) 



-AiWa 2 ( Pl -p 2 )([ Pl xp 2 ] -s)(10) 



Here Pi — p 2 = ±Gj, and s is the effective spin 1/2 
introduced in ([3]). At long wavelengths the leading order 
contribution for the SOC is given by: 



S((2\U\1)) 



= -4iWa 2 {K 1 ■ K 2 ){[K l x K 2 \ ■ s) 
= ij=\v\s z , \ V \ = l^K*W, 



which can be written as: 



SH h 



V3 




Projecting this matrix to the states \a) and \b) defined 
by ([7]) we get: SHi s — > — 2\rj\s z <j y . Thus, the final Hamil- 
tonian, including ([5]), reads: 



H = v (~<7 z q x - a x q y ) - 2\r)\s z a y . 



(11) 
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After a unitary transformation the Hamiltonian can be 
written in its conventional form jij: H — > vcr-q—2\r]\s z (j z 
leading to a SOC gap given by: 



tential to the lattice potential ((4]). We assume: 



2|r?| = 3a 2 K 4 W 



16 
~9~ 



W , 



(12) 



which shows that the SOC gap is a very strong function 
of the d/L ratio. 

SOC matrix elements (|10[) can be easily included in 
the exact diagonalization. Notice that eqs. ([3]) and (fTU)) 
are derived assuming a 2 p A <C 1. Even at small a this 
condition is violated for high momenta states included 
in the exact diagonalization. To avoid this problem we 
account (|10[) only for three lowest quantum states and 
set the SOC matrix element equal to zero for all higher 
states. This procedure gives the correct gap near the 
Dirac points. Notice that, a 2 K 4 = if j£ w 0.037 , for 
d/L = 0.5. The dispersion calculated numerically for this 
value of the SOC and W = 2 is shown in Fig. The 
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FIG. 3: (a) The hole dispersion at W = 2 along the T -> 
/G — > -^"3 — ^ r contour in the BZ. The spin orbit split- 
ting corresponds to d/L = 0.5. (b) The dashed lines show 
the dispersion of s z — ±| edge states. The solid lines show 
boundaries of the 2D continuum. 



calculated spin orbit gap is close to the analytical formula 
(fT2| . By varying the transverse confinement width d, see 
Eq . (fT2"j) . one can continuously change the SOC gap and, 
hence, the electronic properties of devices made in these 
systems. 

In order to study the presence of edge states we have to 
consider a finite sample with a confining potential at the 
edge. Once again, having in mind simplicity, we consider 
an infinite edge potential and by adjusting the position 
of the wall we can create edge states. The mechanism of 
the edge state formation discussed here is qualitatively 
different from that in the graphene [H[ because of the 
nature and strength of the potentials in the two cases: 
graphene is better described by the tight-binding model, 
while semiconductors are better described by the nearly 
free electron approximation. In the tight binding ap- 
proximation edge states exist even without the SOC (say, 
along the zig-zag direction), and the SOC only modifies 
the dispersion of the state [4(. On the other hand, in the 
nearly free electron approximation the SOC is crucial for 
the formation of the edge state. 

Let us consider the effect of a laterally confining po- 



if y > y 
oo if y < yo 



The envelope wave function of the edge state at y > yo 
is given by: 



ip = A 



4>a 



-Ay 



Solving Hip = eip with Hamiltonian (fTTj) one finds: 



e = ±y / r] 2 + v 2 q 2 — v 2 X 2 

= 1 

.vX + n 



(13) 



vqx - e 

The corresponding full coordinate wave function reads: 



ipoce lq * x e- Ay {\a)+i> b \b}} 



|6, = -L ( 2e 



iK-2-r giJf 3 r\ 



iKfr _ e iK 2-r _ e iK 3 r\ 



(14) 



where we have used Eqs. ([6]), (JT)) for \a) and \b). The 
wave function must be zero at y = yo at the position 
of the confining wall. It is obvious that one cannot sat- 
isfy this condition at arbitrary yo. Hence, generically 
the edge state either does not exist or requires a more 
involved calculation to find it. The calculation must in- 
clude admixture of high momentum states to the wave 
function (fT^|) . Fortunately, it is very easy to find the 
state if the wall position is chosen as: r^^J/o = n N , 
where N is integer. At y = yo the basis wave function 
| a) is zero at any x. Therefore, to have ip{x,yo) = we 
need only to impose ipb = 0. Hence, using (|13p we con- 
clude that the edge state exists only at negative ry, i.e at 
s z = +1/2 (see Eq. @T)): 

s z = +1/2, Tj < : X = -rj/v, e=-vq x , 

which is valid near the K Dirac point. Notice that the 
SOC correction (pTJ|) - ([TT]) . is even under the reflection. 
Therefore, at the K' Dirac point the effective Hamilto- 
nian differs from (|11[) only by the replacement v — > — v . 
The same is valid for the edge state solution ([IB"]) . There- 
fore, at K' the edge state exists only at positive rj, i.e at 
s z = -1/2 (see Eq. (IT])): 



-1/2, r) > : A = rj/v , 



vq x 



The dispersion of the edge states with s z — ±4 is shown 
in Fig. (lb). 

The edge states support the spin current at the edge 
of system in the regime of a TI. This is true if the edge is 
parallel to a face of the BZ. A sample must have at least 



■5 



two confining walls. If the walls are created by external 
gates they can be moved. In the regime of a TI one can 
tune positions of the walls in such a way that there will 
be a spin current on one edge and there is no current on 
another edge. Other combinations are also possible 

In conclusion, we have shown that it is possible to 
create robust TI states in ordinary semiconductors with 
strong SOC by quantum confinement and superimposed 
potentials with hexagonal symmetry. We have shown 
that the SOC gaps can be as large as the bandwidth and 
can be controlled by varying the confinement potential, 
the strength and scale of the superimposed potentials. 
These systems present amazing flexibility and can be 
tuned between completely different regimes such as the 
Dirac metal, TIs, and standard band insulators. Thus, 
they present an opportunity to study exotic physics in 
the framework of materials that have been the basis of 
the current semiconductor technology. 
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